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Abstract 

We study direct detection of dark matter in supersymmetric (SUSY) model 
where most of SUSY particles have very high-scale masses beyond the weak scale. 
In the scenario, Wino-like or Higgsino-like neutralino is a good candidate for the 
dark matter in the universe. The neutralino scatters off nuclei by a Higgs-boson 
exchange diagram and also elect roweak loop diagrams. It is found that the elastic- 
scattering cross section with nuclei is enhanced or suppressed due to constructive 
or deconstructive interference among the diagrams. Such a cross section is within 
the reach of future experiment in some parameter region. 



1 Introduction 



Supersymmetric (SUSY) remodelling of the Standard Model (SM) is one of the promis- 
ing candidates for physics beyond the SM. The minimal extension, called the minimal 
supersymmetric Standard Model (MSSM), has been studied enthusiastically in various 
literature. The weak-scale SUSY is, however, severely constrained by the experiments 
at the Large Hadron Collider (LHC). Since no signal of SUSY particles has been discov- 
ered yet, the ATLAS and the CMS Collaborations have imposed stringent limits on their 
masses, especially those of colored particles |1|. The weak-scale SUSY is also challenged 
by the discovery of the SM-like Higgs boson with a mass of about 125 GeV, which is re- 
cently reported by the collaborations [2l . In the MSSM radiative corrections from heavy 
sfermions make the Higgs mass larger |3|. Thus those results from the LHC may indicate 
that the SUSY scale is somewhat higher than the weak scale El. 

Although the high-scale SUSY scenario sounds unnatural in a viewpoint of the hier- 
archy problem, phenomenological aspects of heavy sfermions are quite fascinating |5]jll 
Due to the sufficient radiative corrections, the 125 GeV Higgs boson may be achieved [12 
The SUSY contributions to flavor changing neutral current (FCNC) processes and electric 
dipole moments are suppressed by heavy sfermion masses so that the SUSY flavor and CP 
problems are relaxed 113) . In cosmology, the gravitino problem may be avoided because it 
may be as heavy as sfermions, then the thermal leptogenesis for baryon asymmetry in the 



universe works with high reheating temperature 14 . On top of that, the gauge coupling 



unification is achieved as precisely as that in the MSSM since the sfermions form the 
SU(5) multiplets, and the proton lifetime could be well above the current experimental 
limit (15). These features have stimulated various works |16|. 

The high-scale SUSY scenario does not necessarily mean all SUSY particles in the 
MSSM are heavy. Superpartners of gauge bosons and Higgs bosons may be at the weak 
scale without destroying above features. This is plausible because the lightest particle 
among their mixed states is the lightest SUSY particle (LSP) and it is a good candidate 
for the dark matter (DM) in the universe. Such a candidate is one of so-called Weakly 
Interacting Massive particles (WIMPs). Though sfermions may be beyond the reach of 
the LHC, the LSP DM may be searched in the direct dark-matter detection experiments. 

In this article we consider a scenario where the LSP mass is around the weak scale 
and the other SUSY particles are much heavier, and give a precise calculation of elastic- 
scattering cross section between the LSP DM and nucleon. Wino-like or Higgsino-like 
neutralino is a viable DM candidate since the thermal production in the early hot uni- 
verse gives the observed DM density even in the heavy sfermion scenario; Wino with a 
mass of 2.7-3.0 TeV [17] or Higgsino with a mass of 1 TeV |18|. The neutralino mass less 
than TeV is also possible to explain the the DM density when its non-thermal produc- 



tion is considered [L9, 20 . Also the Wino LSP is a natural consequence of the anomaly 
mediation 121). Therefore we focus on those well-motivated cases. Since there are only 
a few undetermined parameters, the observed value of the Higgs-boson mass at 125 GeV 
allows us to make a robust prediction for the scattering cross section with nucleon. As 
we will see below, both the tree-level |22,23| and the loop-level processes [241426] give rise 
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to sizable contributions to the scattering cross section. 

This paper is organized as follows. In the next section, we explain the scenario of high- 
scale supersymmetry, in which the Wino-like or Higgsino-like neutralino is predicted as the 
DM. In Section 3, the effective Lagrangian for the neutralino- nucleon elastic scattering is 
reviewed. In Section 4, we discuss relevant tree and loop-level contributions to the spin- 
independent (SI) scattering of the Wino-like or Higgsino-like neutralino, and evaluate 
the cross section. Section 5 is devoted to conclusion. In the evaluation of the SI cross 
section in text, we use the results of the lattice QCD simulations for the mass fractions 
of light quarks in nucleon. We also show the results when the mass fractions estimated 
in the chiral perturbation theory are used in Appendix A. In Appendix B, we give the 
spin-dependent (SD) cross section for completeness. 



2 The Scenario 

In this section we briefly describe the scenario which we discuss in this paper. As it is 
mentioned in the Introduction, we consider a SUSY scenario where all SUSY particles 
are well above the weak scale, except for Wino or Higgsino. Such mass spectrum is given 



by a simple SUSY breaking mechanism [5j|6, 11,21|. Assume that there exists a SUSY 
breaking hidden sector containing a SUSY breaking field Z which is charged under some 
symmetry. Then a generic form of Kahaler potential yields masses of Msusy ~ F z /M* 
for all the scalar bosons in the MSSM except the lightest Higgs boson (F z and M* are 
the F-component vacuum expectation value (VEV) of the field Z and the messenger 
scale, respectively). On the other hand, since Z is charged under some symmetry, the 
gaugino and Higgsino mass terms are not given by the Z-field linear terms. Thus they 
do not necessarily have the mass scale Msusy, i.e., they are model-dependent Q Now we 
consider the case of M* = Mpi (Mpi is the reduced Planck scale). In the case, gauginos 
acquire their masses via the anomaly mediation, which are of the order of /167r 2 . 
Here m 3 / 2 = Fz/VSMp\ is the gravitino mass. For Higgsino, on the other hand, so-called 
\i term, /LiH u Hd in the superpotential (H u and are up-type and down-type Higgs 
chiral superfields, respectively), may be absent by a certain symmetry, e.g., the Peccei- 
Quinn symmetry [27] . In such a case, the gaugino-Higgs loops induce the Higgsino mass, 
which is smaller than the gaugino masses by another loop factor. Thus, the Higgsino-like 
neutralino is the LSP. On the contrary it may be as heavy as gravitino in another case. 
When the Kahaler potential has a term, K = tvH u Hd + - • • , the Higgsino mass is provided 
by the supergravity effects and it lies around the gravitino mass scale. In that case, the 
Wino-like neutralino is the LSP. 

In order to consider the Wino-like or Higgsino-like neutralino, we take the Wino and 
Higgsino mass parameters as free parameters. Namely, the gaugino and Higgsino mass 



1 Trilinear soft SUSY breaking terms are also model-dependent. However, they are irrelevant in our 
discussion. 
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terms are given by 

Co = " E \M a KK-^H u H d , (1) 

a=l,2,3 

where H u H d = H+Hj - H*H° d and 

M a = ^m 3/2 . (2) 

Here X a (a = 1,2 and 3) are Bino, Wino and Gluino, respectively, and H u and Ha are 
up-type and down- type Higgsinos, respectively. The coefficients b a denote the one-loop 
beta-functions of the gauge coupling constants g a (a = 1,2 and 3 for U(l)y, SU(2) ^ 
and SU(3)cj respectively) given by (&i, 62,^3) = (33/5, 1, — 3)0 After the electroweak 
symmetry breaking, Bino (5), Wino (VF ), and neutral Higgsinos mix with each other. 
The mass eigenstates, called neutralinos, are obtained as Xi = J2j %ijXp where Xi = 
B,W°,H® and H% for i — 1,2,3 and 4, respectively. Xi is the lightest neutralino, and 
from now on, we omit the subscript of Xi f° r simplicity. In the following calculation we 
take M 2 as a free parameter instead of m 3 / 2 , and consider the case where the lightest 
neutralino explains the current relic density of DM. 

When sfermions are very heavy, the 125 GeV SM-like Higgs boson is achieved with 
tan/3 r>u 1-5 |12|. Here tan/3 is the ratio of the VEVs of up- and down-type Higgs fields. 
In this article, however, we also consider larger tan/3 to demonstrate the cross section 
for general case, assuming appropriate sfermion masses to make the Higgs boson mass 
125 GeV. 



3 Neutralino-Nucleon Scattering Cross Section 

Here we give formulae for the calculation of the scattering cross section of the neutralino 
with nucleon [25,28,29). It is calculated from the effective Lagrangian for scattering of 
the neutralino with quarks and gluon in the limit of low relative velocity, which is given 

by 

£eff = d q xV75X° Tiiil^ + f q m q X°X° QQ 

(1) (2) 

+ \ X°id^X° 0% + |L_ x°{id^)x° 0%, 

(1) (2) 

+ fa X°X°G%G a ^ + ^ x°id^X° 0% + g_ 0% , (3) 

where M and m q are the masses of the neutralino and quarks, respectively. Sum over 
quark flavors q = u,d,s for the first and second terms and q = d, s, c, 6 for the third 



2 Threshold correction may change the ratio M1/M2 from the above relation. However, this does not 
affect our numerical result significantly. 
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and fourth terms is implicit. The field strength tensor of the gluon field is denoted by 
G a . The last two lines include the quark and gluon twist-2 operators, 0^ u and 0^ 
respectively, which are defined as, 

0% = G<?G% + (g»„/4)G a a pG^ , (4) 

with = <9 M — ig^A^Ta the covariant derivative (T a is the generator of SU(3)c, and A a 
is the gluon field). In order to remove the redundant terms, we use the integration by 
parts and the equation of motion for the operators. The first term in Eq. ^ yields the 
SD interaction, while the other terms generate the SI interactions. 

In order to compute the x°-nucleon cross section from the effective Lagrangian, we 
need to evaluate the nucleon matrix elements of the quark and gluon operators. The 
nucleon matrix elements of the scalar-type light-quark operators, i.e., m q qq (q = d, s), 
are parametrized as 

(N\m q qq\N) = m N f^ , (5) 

where is the nucleon mass and \N) denotes the one-particle state of the nucleon 
(TV = p, n). For the heavy quarks and gluon, on the other hand, their matrix elements 
are obtained by using the trace anomaly of the energy-momentum tensor in QCD: 

— cy 
(N\m Q QQ\N) = -r^c Q {N\Gl u G a ^\N) , 

rnsf^ = ~(N\G^G^\N) , (6) 

with f^} = 1 — Y^qfrq^ an d a s = dl/^- The long-distance QCD correction cq in the 
above expression is evaluated in Ref. [30] as cq = 1 + lla s (mg)/47r, and we take their 
numerical values as c c — 1.32, = 1.19 and c t — 1 in this paper. As can be seen from 
Eq. ^ , the scalar-type heavy quark operators contribute to the nucleon matrix elements 
only through the loop-induced gluon operator. 

The nucleon matrix elements of the twist-2 operators are evaluated with the parton 
distribution functions (PDFs): 

(N(k)\0%\N(k)) = —(k,K - m 2 N g^/4) (q N (2) + q N (2)) , (7) 
m N 

{N(k)\0%\N(k)) = —{k,k v - m 2 N g^/A) G N (2) , (8) 
m N 

where ##(2), ^(2) and Gn{2) are the second moments of PDFs of quark, anti-quark and 
gluon, respectively, which are given by 



q N (2) + q N (2) = dx x [q N (x) + q N (x)] 

Jo 
pi 



G N {2) = I dx x gjsr(x) . (9) 



o 
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Here, we use the PDFs at the scale of \i — mz {rnz is the Z boson mass), since, as will 
be described later, the terms with quark twist- 2 operators in Eq. ^ are induced by the 
one-loop diagrams in which the loop momentum around the weak boson mass scale yields 
dominant contribution. 

Finally, the SI effective coupling is obtained as 

A = fW f -U - (a„(0\ + a„(0\\ (nW + - — fW f„ + — H^O.} ( + 



C/ 9 + \ (qs{2) + q N (2)) (gW + <£>) - ^f®fo + \g n {2) ($ + «£>) . (10) 



Here the sum of quark flavors is implicit as in Eq. ([3]). Note the factor l/a s in front 
of fc in Eq. (10). It makes the gluon contribution sizable, although the interactions 
of the neutralino with gluon are induced by higher-loop processes than those with light 
quarks [25]. On the other hand, the contributions of the twist-2 operators of gluon are 
sub-dominantj^] as g<p and g^? are suppressed by the strong coupling constant a s . Thus, 
we ignore them in this paper. 

The effective axial vector coupling, which is relevant for the SD cross section, is readily 
written as 

a N = d q Aq N , (11) 

with 

(N\q^q\N) = 2 SfI Aq N . (12) 

Here is the spin of the nucleon and quark flavor sum is taken for q = u,d, s. By using 
the effective couplings obtained above, we obtain the cross section of the neutralino with 
nucleon: 

a N = -m 2 R [\f N \ 2 + 3\a N \ 2 ] , (13) 

7T 

where ttir = Mm^j[M + mjsr) is the reduced mass of neutralino- nucleon system. 

Before concluding this section, we refer to the numerical values for the parameters 
that we use in this paper. The mass fractions of light quarks, defined in Eq. (5) 



are to be extracted from the results of the lattice QCD simulations 32,33). The mass 
fractions of light quarks are evaluated with independent methods and also by independent 
groups so that the results derived with the lattice QCD simulations have become more 
reliable. The mass fractions evaluated in the chiral perturbation theory (ChPT) have 
larger uncertainties than those from the lattice QCD. The SI cross section evaluated by 
the use of the mass fractions from the ChPT, which predicts larger is shown in 



3 Evaluating the nucleon matrix elements of twist-2 operators at \i = mz makes the perturbative 
expansion with respect to a s reliable. Instead, if one would like to estimate the matrix elements at 
fx = 1 GeV, one also needs to take into account the operator-mixing effects due to the QCD radiative 
corrections, and therefore, to include the gluon twist-2 operator, as in Ref. [3l]. These two approaches 
are equivalent since the sum of the terms with twist-2 operators is scale-independent once they are 
multiplied by their coefficients. Generally speaking, however, the former approach makes the calculation 
robust thanks to the perturbativity of the QCD coupling. 
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Mass fraction 


(l 


3roton) 


Ap) 

JTu 

Ap) 

JTd 

Ap) 

JTs 


0.019(5) 
0.027(6) 
0.009(22) 


(neutron) 


An) 
JTu 
An) 
JTd 
An) 
JTs 


0.013(3) 
0.040(9) 
0.009(22) 



Second moment at ji = m z 


u{2) 


0.22 


«(2) 


0.034 


d(2) 


0.11 


d(2) 


0.036 


s{2) 


0.026 


5(2) 


0.026 


c(2) 


0.019 


c(2) 


0.019 


6(2) 


0.012 


6(2) 


0.012 



Spin fraction 


Au p 


0.77 


Ad p 


-0.49 


As p 


-0.15 



Table 1: Parameters for quark and gluon matrix elements. Errors are shown only for the 
mass fractions, which are used for comparison with the cross section evaluated with the 
mass fraction from the ChPT. 



Appendix A for comparison. The second moments of the PDFs of quarks and anti-quarks 

The spin fractions, Aq N 



m 



are calculated using the CTEQ parton distribution 34 
Eq. (12) are obtained from Ref. [35]. In Table [l] we list the numerical values for the mass 
fractions of both proton and neutron as well as the second moments of the PDFs and the 
spin fractions for proton. The second moments and the spin fractions for neutron are to 
be obtained by exchanging the values of up quark for those of down quark. 



4 Results 



In this section we calculate the SI cross section for the Wino-like or Higgsino-like neu- 
tralino in the high-scale SUSY scenario. The values of the Higgsino and Wino mass 
parameters are model-dependent. Therefore, we regard both M 2 and ji as free parameters 
in the following analysis while we take M 2 positive^] 



The tree-level x X -Higgs interaction yields scalar-type effective operators, x^X^QQ 



and X^y^G^G^ (Left in Fig. 1). Let us denote these contributions by and /g, 
respectively. By adding them to the loop-level contribution^] J^wimp an j jewimp w j 1 i c ] 1 



4 We assume the parameters to be real in this article. Possible phases of the parameters might affect 
the "Higgs" contribution, which is defined later, to the SI effective coupling. 

5 Those contributions are evaluated in a pure Wino or Higgsino limit. When — /i, one needs 
to take the mixing among them into account and modify the formulae in Ref. p6] appropriately. In 
the present situation, however, the tree-level contributions dominate the loop-loop level ones. Thus, the 
modification in the loop-level effects has no significance on the resultant scattering cross section. There 
exists a case where the tree-level contribution is still subdominant even when — \i. As we will see 
later, however, the lightest neutralino is almost pure gauge eigenstate in such a case. Thus the results in 



Ref. 26 are applicable. 
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Higgs 



r\/\/\/\/\j 

z 



axial-tree 



Figure 1: Diagrams via tree-level x°-x°-Higgs/Z interaction in elastic x°-nucleon scat- 



tering. "Higgs" contribution and "axial-tree" contribution are defined in Eqs. (21) and 

pi. 



are induced via the W/Z boson loop diagrams |26|, we obtain 



u 

fa 



jEWIMP 

rH , fEWIMP 
JG "T" JG ? 



with 



rH 
Jq 

rH 
JG 



gls h 

Am w m z h 
a 



2 ' 



12tt 



H 



(14) 
(15) 



(16) 
(17) 



where mw and are the masses for W boson and Higgs boson, respectively. The 
coupling of the neutralino with the Higgs boson is denoted by s h in the above expression, 
which is given as 



s h = (Zyi — Zn tan 9 W ) (Z\% cos — Z\± sin 0) . 



(18) 



Here Ow is the weak mixing angle and we take the decoupling limit since the heavier Higgs 
bosons have masses much larger than the weak scale. In addition, when M 2 , > m\y, 
the coupling s h is approximated as 



s h ~ 



in the Wino-like neutralino case and 



(M 2 + fi sin 2/3), 



s h ~ 



m w 



Mo 



m w tan 2 9 W 
Mi - \n\ 



(1 ± sin 2/3), 



(19) 



(20) 



in the Higgsino-like neutralino case. Here the plus (minus) sign in front of sin 2/3 is for 
/x>O(/x<O)0 



3 There is an sign error in Eq. (24) of Ref. 



the correct expression is d=i 



36 



raw 



M 2 -|/x| 



mw tan Ow 
Mi-H 



for /i < case. In addition, for heavy Higgs coupling, 
cos 2j3. 
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r\/\/\/\/\j 

w/z 

q 
q 

w/z 



scalar, twist-2 
axial-loop 



scalar 



HQ 




sduon 



duon 



Figure 2: Diagrams which are induced by electroweak interaction in elastic x°-nucleon 
scattering, "scalar" , "twist-2" , "gluon" and "axial-loop" correspond to each term in the 



effective couplings. Their definitions are given in Eqs. (21) and (24). A complete set of 
diagrams is given in Ref. |26|. 



As it is seen in Eqs. (19) and (20), in the case where one mass parameter is much larger 
than the other (i.e., M 2 <C or <C M 2 ), the lightest neutralino becomes almost pure 
Wino or Higgsino state. Then the tree-level x°-x°-Higgs interaction, as well as 
interaction which is relevant for the SD scattering, is suppressed. Thus the loop-level pro- 
cesses become important. The loop-level effective couplings are calculated in the previous 
work 1 26 1, where the elastic scattering cross section for generic electroweak-interacting DM 
particles (i.e., n-tuplet of SU(2) L with hypercharge Y of U(l)y) is evaluated. Pure Wino 
corresponds to n = 3 and Y — 0, while pure Higgsino corresponds to n = 2 with Y — 1/2. 
The previous results have revealed that the loop-level contributions are sizable when the 
DM-particle mass is much larger than those of weak bosons. (See also Ref. [36].) Further, 
it has been found that the SI cross section tends to be suppressed with the 125 GeV Higgs 
boson mass due to an accidental cancellation. (See Fig. 5 of Ref. [26].) These observations 
indicate that both the tree-level and the loop-level contributions are significant in a wide 
range of parameter space. Taking the above discussion into account, we calculate the 
scattering cross section of the neutralino with nucleon including all the possibly dominant 
contributions. 

For later discussion, we refer to each term in Eq. ^ as, 

fN AN) f EWIMP , 3 / (0 \ i - / 9 \\ (Jl) | „(2)\ 87F AN) f EWIMP 

— ~ JTq Jq +4 KQnW + Qn{^)) [g y q + 9\ ) ~ J TG J G 

8tt 
9a 



m N 



s 



I / AN) rH _ AJN) r 

^ I J Tq Jq q JTG J G 



[(scalar) + (twist-2) + (gluon) + (Higgs)] /m N . (21) 
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Here "scalar" , "twist-2" and "gluon" contributions are from diagrams in Fig. [2| while we 
define "Higgs" contribution, which contains both and as shown in Fig. [lj 

Before turning to numerical calculations, we briefly discuss the gaugino-siermion- 
fermion couplings, which we call the gaugino couplings hereafter, in the high-scale SUSY 
scenario. The gaugino couplings are equal to the gauge couplings at the energy scale 
larger than Msusy- With the scalar particles decoupled at Msusy, however, the effective 
theory below the scale is not supersymmetric anymore; thus the gaugino couplings might 
in general deviate from the relations. This deviation affects the neutralino mass matrix, 
leading to corrections to Zij. Using the renormalization group equations for the gaugino 
couplings in the split SUSY scenario given in Ref . [7] , we explicitly calculate the running 
of the couplings and find that the deviation of gaugino couplings from the corresponding 
gauge couplings is less than a few %; e.g., the U(l)y gaugino coupling decreases from 
the supersymmetric one by around 7%, while the SU(2) L gaugino coupling increases by 
about 1%, when Msusy is 10 3 TeV and tan/3 = 1, which gives the Higgs mass of around 
125 GeV. 

Now we are ready to give numerical results of the scattering cross section. Fig. [3] shows 
the results of the SI scattering for the Wino-like neutralino. Here we give the cross section 
of the neutralino with proton. In the plots we take M 2 = 3 TeV and = 125 GeV and 
\i < (/z > 0) for the top (bottom) panel, tan/3 = 1.1, 2 and 50 are taken in left, middle 
and right panels, respectively^ In the plot of the SI cross section, green dashed lines 
indicate the SI cross section with only the Higgs contribution taken into consideration, 
while purple solid lines show the result with all the leading contributions included. Shaded 
regions imply error coming from the mass fractions. It is found that the loop contribution 
is important in wide range of parameter space. Let us see \i < case first. When tan /3 < 2, 



the Higgs contribution scales as oc s h ~ m^/(M 2 + M) from Eq. (19), which does not 
grow as | ji \ — M 2 gets smaller. As a consequence, the loop contribution is comparable or 
larger than the Higgs contribution, depending on — M 2 . This behavior is seen in the 
plot of f p in Fig. [ij Here we give the plot of each contribution defined in Eq. (21). In 
the plot result for tan /3 = 50 is given for scalar, twist-2 and gluon contributions (though 
they are insensitive to tan/3), and results for tan/3 = 1.1, 2, 5 and 50 are given for the 
Higgs contribution from top to bottom. Since the Higgs contribution is always positive for 
tan/3 ~ 1.1, it interferes constructively to the other contributions. When tan/3 3> 1, on 
the other hand, situation gets changed. In this case the Higgs contribution is sensitive to 
M 2 and — M 2 for given tan /3. This fact can be seen from Eq. (19). Now s h is given as 



s h ~ mw M2 M2_ J? n 13 • Thus the Higgs contribution is negative when < M 2 tan /3/2 and 



flips its sign in larger \ji\ region. This causes an cancellation; as — M 2 gets larger, the 
absolute value of the tree-level Higgs exchanging contribution drops, and all the negative 
contributions cancel the positive twist-2 contribution around |/i| — M 2 ~ a few hundred 



7 Note that when tan f3 = 1 the tree-level axial coupling, as well as the tree-level Higgs coupling for 
the Higgsino-like DM for negative /i, vanishes exactly. However, it is not the realistic case. In fact the 
gaugino coupling in the neutralino mass matrix receives a correction from renormalization-group effects 
from high scale, which we discussed above, and as a consequence the tree-level couplings do not vanish. 
For the purpose of studying tree- and loop- level contributions in general, we simply avoid tan/3 = 1. 
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GeV to a few TeV, depending on tan/3. We have checked that the result is almost the 
same when tan /3 > 5. 

When /i > 0, on the contrary, the Higgs contribution is always negative. The right 
panel in Fig. [4] shows it for tan/3 = 1.1, 2, 5 and 50 from bottom to top. Thus resultant 
cross section is similar to those in the case of large tan/3 and \i < 0. In both \i < 
and ji > cases, the Higgs contribution becomes irrelevant and the loop contribution 
dominates the cross section when \fi\ — M 2 > 10 TeV. Then the SI cross section lies 



around a value of ~ 10 -47 cm 2 , which is consistent with the results in Refs. [24-26 . 

We have checked the SI cross section is almost independent of the neutralino mass 
except for \i < and low tan/3. We also give the result for M 2 = 200 GeV in Figs. [5] 
and [6j Here we take the other parameters the same as those in Fig. [3j In this case, the 
Higgs contribution becomes larger, leading to a bit enhanced SI cross section. However, 
with relatively large tan/3, a cancellation happens then the cross section behaves similar 
to the previous results, as it is seen in the figure. 

Next let us discuss the Higgsino-like neutralino case. The results are shown in Figs. [7] 
and [8j Here we take |/i| = 1 TeV and rrih = 125 GeV. The upper and lower panels 
correspond to /x < and \i > cases, respectively, tan/3 is taken as similarly to that in 
Figs. [3]and[4j In the fi < case, the Higgs contribution is suppressed by (1 — sin 2/3) for 
tan/3 ~ 1, then the loop contributions become dominant. This is clearly seen in Fig. [5J 
The cross section is around 10 -49 cm 2 in the region M 2 — |/i| > 500 GeV. When tan/3 
is larger, the Higgs contribution scales as s h ~ —mwj[M2 — M) and becomes dominant 
in the effective coupling in the region M 2 — |/i| < a few TeV to 10 TeV, depending on 
tan /3. In this case a cancellation occurs around M 2 — ~- a few dozens of TeV, and the 
SI cross section is about 10 -49 cm 2 for larger values of M 2 — Similar cancellation is 
observed for the \i > case. In the \i > case, the Higgs contribution is not suppressed 
around tan /3 ~ 1 in contrast to \i < case. That is why a significant cancellation always 
happens. 

Here we briefly comment on contributions by the heavy Higgs boson, which we did not 
take into account. The heavy Higgs-x°-X° coupling s H is given by s H ~ — A ^^ 2 cos 2/3 for 

the Wino-like neutralino and s H ~ ±\ + m ^f_ n | ^ l cos 2/3 for the Higgsino-like 

neutralino. Here the overall positive and negative signs correspond to the \i > and 
H < cases, respectively. As it is seen, contributions from the heavy Higgs boson are 
suppressed when tan/3 1. Even when tan/3 > 1, it is suppressed by the heavy Higgs 
mass. 

Finally, we give contour plots of the SI cross section on the — M 2 plane in Fig. [9j 
In the plot we take = 125 GeV. Upper and lower panels are for the cases of fi < 
and ji > 0, respectively, and tan/3 is taken as 1.1, 2 and 50 from left to right in each 
panel. Purple solid lines are contours of the SI cross section of full calculation and green 
dashed lines show the ones given by the Higgs contribution only. The contours of the cross 
section smaller than 10 -49 cm 2 are not shown here. In the figure, dark-shaded regions 



are excluded by the XENON100 experiment 37 . The blue dot-dashed lines and light 



shaded regions correspond to prospected reaches of future experiments. To evaluate those 
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sensitivity limit, we use a value of 10 -47 cm 2 and 10 -48 cm 2 at a DM-particle mass of 
60 GeV and rescale it with respect to the DM-particle mass. The former value is based 
on a discovery sensitivity in a ton- year experiment^] while the latter comes from the fact 
that a sensitivity for the cross section of less than 10 -48 cm 2 is difficult to be achieved 
due to atmospheric neutrino background [38] . 

In the Wino-like neutralino region, it is seen that the full calculation deviates from 
the one given by the Higgs contribution significantly. This is due to the suppression of 
the Higgs contribution (especially for the /i < and tan /3 ~ 1 case) or the cancellation 
in the effective coupling. In addition, as we discussed previously, the cross section can be 
enhanced due to constructive interference between the Higgs and twist-2 contributions in 
the \i < 0, M 2 ^ TeV and low tan/3 region (See left and middle in upper panel). Thus 
such a region can be probed in the future experiment even when is as large as a dozens 
of TeV. If much better sensitivity was accomplished, larger could be also studied. 

For the Higgsino-like case, the Higgs contribution almost determines the cross section 
in the region where future experiments may reach. The loop effect becomes important 
when the the Higgs contribution is suppressed (for /i < and tan f3 ~ 1 case) or in the 
region where M 2 is above several TeV. For the \i < and tan /5 c± 1 case, the cross section 
is around 10 -49 cm 2 thus it is far below the sensitivity of future experiments. 



5 Conclusion 

The high-scale supersymmetry, in which SUSY particles are much heavier than the weak 
scale except for gauginos and/or Higgsino, is favored from viewpoints of the discovered 
125 GeV Higgs boson, null results in the SUSY particle searches at the LHC, and the 
SUSY FCNC and CP problems. In this scenario, the Wino-like or Higgsino-like neutralino 
is a good candidate for the dark matter in the universe. While Wino with a mass of 2.7- 
3.0 TeV or Higgsino with a mass of 1 TeV is predicted in the thermal relic scenario, the 
non-thermal production may explain the observed dark matter abundance in even lighter 
mass. In this scenario, the elastic scattering of the neutralino with nucleon, relevant to 
the direct dark matter search experiments, is induced by tree-level Higgs boson exchange 
diagrams and also loop diagrams due to the electroweak interaction. 

In this article, we evaluate the SI cross section of the Wino-like or Higgsino-like neu- 
tralino including contribution from the loop diagrams of the weak bosons. Since the loop 
diagrams are not suppressed by the neutralino mass, they may be comparable to or even 
dominate over the Higgs-exchange contribution, especially for the Wino-like neutralino. 
As a result, the SI cross section is sensitive to the sign of \± and tan/5, in addition to 
absolute values of the Higgsino and Wino mass parameters since the diagrams are con- 
structively or destructively interfered with each other. 

Due to atmospheric neutrino background, it is difficult to discover the DM in the direct 
detection experiments when the SI cross section is smaller than 10 -48 cm 2 at a DM-particle 



For reference, 8.6 x 10 cm at a DM-particle mass of ~ 60 GeV is the sensitivity of 90% C.L. 



discovery at a ton- year Xenon target experiment 



38 
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mass around 60 GeV. We found that the prediction for the SI cross section is larger than 
the limit in broad parameter region. (See Fig. |9j) The large-scale experiments for the 
direct DM detection are hopeful. 
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A Spin-independent Cross Section Evaluated with 
Mass Fraction from Chiral Perturbation 

In this paper we have used the input parameters extracted from the lattice QCD sim- 
ulations for the mass fractions f^. As a result, the error of the calculation is small, 
as we have seen above. However, another result is also reported for the mass fractions 
based on the chiral perturbation theory. In this case the mass fractions for proton are 



given as = 0.024(4), f$ = 0.041(6) and fj$ = 0.40(14) Large discrepancy 

is seen for as well as larger error. To see the impact of the mass fractions on the 
SI cross section, we plot the results using the mass fractions extracted from the ChPT 



(and the other parameters are unchanged) in Figs. 10 and 11 for the Wino-like case, and 



Figs. [12] and [13] for the Higgsino-like case. In both cases we find that the theoretical 
error of the cross section is much larger than those presented in Figs. [3]and[7| Let us see 
the Wino-like neutralino case, for example. The SI cross section has error of an order of 
magnitude when — M 2 > 10 TeV (a few TeV) for tan/3 = 1.1 (50), and what is worse, 
the lower value is undetermined for the larger values of In the Higgsino-like case, it is 
seen the error is much larger than those in the result which is based on the lattice QCD 
simulation. Therefore, we conclude that in using the input of the mass fractions based on 
the ChPT the SI cross section can not be predicted due to the large uncertainty. 



B Spin-dependent Cross Section 

For completeness, we show the results for the SD cross section. In the case of the SD 
scattering, the tree-level axial vector coupling is induced through the Z boson exchange. 
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With the loop-level contribution combined, the axial vector coupling is given as 



rf, = 4 ree + ^ WiMF . (22) 
Here rf^ WIMP is taken from Eq. (4.3) in Ref. [26], and the tree-level contribution is 

dT = ^r(\Ziz\ 2 -\Zu\ 2 )Tl (23) 



with the weak isospin of light quarks. Then as in /jy, we call each term in Eq. (11), 



a N = d! xee Aq N + d EmMP Aq N 

= (axial-tree) + (axial- loop). (24) 

The first term is derived from the W/Z box diagrams shown in Fig. [2j while the second 
term is given in the tree- level Z exchange (the right diagram in Fig. 

The SD cross section for the Wino-like and Higgsino-like cases are presented in Fig. [14 , 
while the effective axial coupling is given in Fig. [15] In the plots we take the same values 
for the SUSY parameters as those in Fig. [3] and Fig. [7] for the Wino-like and Higgsino-like 
cases, respectively. As is obvious from Eq. ( |23| ), the result is independent of the sign of 
[i. In both cases, tree- level and loop- level couplings are constructive. While the tree-level 
contribution highly depends on tan (3 when tan f3 < 10, it turns out to be insensitive to 
tan f3 otherwise. The SD cross section obtained is are so small that there is little hope to 
detect DM via the SD interactions in the future experiments. 
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Wino-like DM (M 2 = 3 TeV, n<0, lattice) 
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Figure 3: SI cross section of Wino-like neutralino with proton. We take M 2 = 3 TeV, 
rrih = 125 GeV and \i < (top) and \i > (bottom), tan f3 is taken as 1.1 (left), 2 (middle) 
and 50 (right) in each panel. Purple solid lines show the result from full calculation, while 
the result only using the Higgs contribution is in green dashed lines (color online). Shaded 
regions show error from the mass fractions of light quarks evaluated with the lattice QCD 
simulation. 



16 




Figure 4: Each contribution in the effective coupling for Wino-like neutralino. ji < 
(left) and \± > (right) and the other parameters are the same as those in Fig. [3j In the 
panel "scalar" (orange dotted), a twist-2" (blue dashed), "gluon" (purple dot-dashed) and 
"Higgs" (green solid) contributions defined in Eq. ( [21] ) are given. For scalar, twist-2 and 
gluon contributions we take tan /3 = 50, while Higgs contribution is given for tan /3 — 1.1, 
2, 5 and 50 from top to bottom (from bottom to top) in the left (right) panel. 
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Wino-like DM (M 2 = 200 GeV, \i<0, lattice) 
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Figure 5: Similar plots to those in Fig. [3] except for taking M 2 = 200 GeV. 
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Wino-like DM (M 2 = 200 GeV, lattice) 
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Figure 6: Similar plots to those in Fig. [4] except for taking M2 = 200 GeV. 
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Figure 7: SI cross section of Higgsino-like neutralino with proton. We take \n\ = 1 TeV, 
rrih = 125 GeV and /i < (top) and /i > (bottom), tan/3 is taken as 1.1 (left), 2 
(middle) and 50 (right) in each panel. Line contents are the same as those in Fig. |3| 
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Figure 8: Each contribution in the effective coupling for Higgsino-like neutralino. \i < 
(left) and \i > (right) and the other parameters as well as line contents are similar to 
those in Fig. [7j 
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Figure 9: Contour of the SI cross section in cm 2 unit. Upper panels are in the case 
where fi < and tan/3 = 1.1 (left), 2 (middle), and 10 (right) are taken, respectively. 
In the lower panels \i is set to be positive. We take rrih = 125 GeV. Results from full 
calculation and the only Higgs contribution are given in purple solid and green dashed 
lines, respectively. Lines are shown for the cross section larger than 10 -48 cm 2 . Here we 
also show the exclusion regions by XENON100 |37| in dark shade. Blue dot-dashed lines 
and light-shaded regions correspond to future prospects of experiments with sensitivities 
to cross section 10 -47 cm 2 and smaller than 10 -48 cm 2 at a DM-particle mass of ~ 60 GeV, 
respectively. 
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Wino-like DM (M 2 = 3 TeV, \i<0, ChPT) 
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Figure 10: Similar plots to those in Fig. [i] except that the mass fractions f^p from the 
ChPT are used. 
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Figure 11: Similar plots to those in Fig. [i] except that the mass fractions f^p from the 
ChPT are used. 
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Figure 12: Similar plots to those in Fig. |7| except that the mass fractions fj, q from the 
ChPT are used. 
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Higgsino-like DM (M = 1 TeV, ChPT) 
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Figure 13: Similar plots to those in Fig. ^ except that the mass fractions from the 
ChPT are used. 
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Wino-like DM (M 2 = 3 TeV) 
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Higgsino-like DM (l|il = 1 TeV) 
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Figure 14: SD cross section for Wino-like neutralino (top) and Higgsino-like neutralino 
(bottom) where the input parameters are set to be the same as those in Figs. [3] and [7| 
respectively. Results are shown for one-loop level and tree-level in purple solid and green 
dashed lines, respectively. 
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Wino-like (M 2 = 3 TeV) Higgsino-like (I |x| = 1 TeV) 
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Figure 15: Each contribution in the effective axial coupling for Wino-like (left) and 
Higgsino-like neutralinos (right). In the figure tree-level and one-loop level contribu- 
tions, defined in Eq. ([24]), are shown in green dashed and orange solid lines, respectively. 
For tree-level contribution tan/3 = 1.1 2, 5 and 50 are taken from bottom to top. 
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